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ABSTRACT 

Results are presented of a numerical calculation of the tunneling gap for a domain wall 
moving in the double well potential of a pair of voids in a magnetic insulator. Both 
symmetric and asymmetric double well potentials are considered. It is found that, even 
in the absence of dissipation, the prospect for observing quantum coherence on a meso- or 
macro-scopic scale appears unlikely. 
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There has been a great deal of interest recently in the prospect that magnetic systems 
might provide a new setting in which to observe a macroscopic degree of freedom behaving 
quantum mechanically [|I|]. To date, the magnetic systems considered are: (i) magnetic 
grains 0, f|, |(J ; and (ii) solitons in magnetic systems (particularly 180°-Bloch walls) 
[1], [7[ ^ 0, |TTJ. Attention has focused on macroscopic quantum tunneling since the 



conditions necessary for its observation appear the most favorable . Arguably, of more 
fundamental interest is macroscopic quantum coherence (MQC) because of its connections 
with quantum measurement theory |13]. In MQC the macroscopic object tunnels periodi- 



cally through the central barrier of a double well potential (DWP). This effect is a direct 
consequence of the quantum state of the object being in a coherent linear superposition 
of macroscopically distinguishable states, and it is this aspect of MQC that connects it to 



the Schrodinger cat paradox [13] and to quantum measurement theory. 

In this paper we will examine quantum coherence (QC) for a 180°-Bloch wall (we sup- 
press the "180°" below) moving in the double well potential due to a pair of voids present 
in a uniaxial magnetic insulator with quality factor Q = K/2ttM 2 ^> 1 (K = magnetic 



anisotropy constant, M = spontaneous magnetization) at T = ||15| . Under these con 



ditions: (i) demagnetization effects are small so that the voids do not alter significantly 



the Bloch wall configuration of the magnetization jig] ; and (ii) most dissipative effects 
are small ||. Thus our system, under the conditions assumed, is expected to have weak 
dissipation — a situation we approximate by ignoring dissipation (though see below). We 
calculate the ground (first excited) state energy E$ (Ei) numerically (from which the tun- 
neling gap is A = E\ — Eq) for: (i) identical voids leading to a symmetric double well 
potential for various wall sizes N (the number of spins in the wall) and void separations 



L; and (ii) non-identical voids leading to an asymmetric double well potential for varying 
degrees of asymmetry (for a particular choice of N and L). We find that observation of 
quantum coherence on either a macro- or meso-scopic scale appears unlikely. For macro- 
scopic QC (N > 10 4 ), weak stray magnetic fields introduce a bias into the gap which masks 
A except for void separations very close to the value at which the central barrier disap- 
pears. At these separations, the tunneling gap varies on a length scale that is less than the 
coarse graining length scale. One would expect that the experimentally relevant gap would 
be a coarse grained average which is seen to be less than the bias and so unobservable. For 
mesoscopic QC (N ~ 10 2 — 10 3 ), although the bias introduced by a stray magnetic field 
is quite small, the slightest bit of asymmetry in the voids is sufficient to pin the wall to 
the larger of the two voids. The difficulty here is that an adequately large tunneling gap 
requires a very small central barrier which is destroyed by the slightest difference in the 
voids. For spherical voids, we find that pinning of the wall can be avoided only when the 
difference in the radii is much less than the coarse graining length scale. Averaging the 
effects of asymmetry over the coarse graining length scale leads one to conclude that a real 
mesoscopic wall will likely be pinned by asymmetry. Thus even in the dissipationless ap- 
proximation, one expects that observation of quantum coherence on any scale larger than 
microscopic appears unlikely due to the severe tolerances imposed on the experimental 
situation. 

The system of interest is a magnetic insulator which is a lattice of spins (with lattice 
constant lo) coupled to each other via the exchange and dipole-dipole interactions; and to 
the underlying lattice via the anisotropy interaction which is assumed to be uniaxial with 
easy axis along z. For the length scales of interest to us the lattice system can be coarse 
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grained so that the magnetic state of the system is described by a magnetization M(x, t) 
defined on a 3-D spatial continuum. The total static energy in the absence of an external 
magnetic field is the sum of the exchange, anisotropy and demagnetization energies. A 
stationary Bloch wall is a soliton configuration of the magnetization M(x) with vanishing 
demagnetization energy, subject to the boundary condition M(y — > ±00) = ±Mz. The 



spatial variation of the wall is localized to a planar region of thickness A = JJ/K (J = 
exchange stiffness constant) which is assumed to be parallel to the xz-plane. The wall 
coordinate q specifies the distance from the origin to a reference point on the wall (viz. 
Xwaii = qy) PI - 

Voids in the magnetic insulator act as pinning sites for the wall. For materials with 
Q > 1, the attraction is due primarily to a reduction in the exchange and anisotropy 
energies that occurs when the wall sits on the void. For a void of length scale R satisfying 
io < < A (io = lattice constant = 5 A), located at the origin, the pinning potential seen 
by a flat Bloch wall is U(q) = — U sech 2 (g/A) |T7|], where U = 2KV d and Vd is the void 



volume. In our calculation, A = 1000A (50 A) for walls with N > 10 4 (300 < N < 3000) 
and U = O.leV (1.0 x l(T 5 eV). For LaGaYIG, Q = 25.2 and K ~ 2000 ergs/cm 3 , so 



x± = ±Qoy (L = 2Q ) with volumes V+ = aV~ (a > 1). They produce the double well 
pinning potential 



When a = 1 we obtain a symmetric double well potential (SDWP); otherwise, the wall 
sees an asymmetric double well potential (AsDWP). For an energy E corresponding to QC 
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that for a spherical void R ~ 200A (IOA). We consider two spherical voids located at 




there will be 4 turning points T\ < T 2 < T3 < T4. We refer to the region q < T\ as the "left 
barrier" ; the region T\ < q < T 2 as the "left well" ; the region T 2 < q < T 3 as the "central 
barrier"; the region T 3 < q < T 4 as the "right well"; and the region T 4 < q as the "right 
barrier". Varying the void separation L (viz. Qq) varies the depth of the wells and the 
height and width of the central barrier. The experimental situation envisioned is either: 
(1) a thin film or; (2) a very narrow wire of the magnetic insulator in which only one Bloch 
wall is present ||. In the thin film case, only a small region (of the wall) of cross-sectional 
area A w is involved in tunneling between the pinning sites ||. In our analysis below, we 
treat the wall as if it were flat, whereas, for the thin film scenario, it will in fact be curved. 
Curvature effects will be discussed below (see also fL6f). In the case of the very thin wire, 
curvature effects are not expected to be important because of the large energy required to 
bend the wall on a length scale of order the cross-sectional dimension of the wire. Flat 



moving walls have a kinetic energy M^q /2 JT5[. Here Md = A w /27ty X is the Doring 



mass; and 7 is the gyromagnetic ratio. The wall Hamiltonian is H = p 2 /2M^ + U(q). 
Introducing the dimensionless length x = q/X and the energy scale S = h 2 /2M D X 2 , we can 
write the time independent Schrodinger equation in the dimensionless form 



d 2 

— j + Uq I sech 2 (a; + xq) + a sech 2 (x — xq) j + S 



V> = 0. 



Here Uq = Uq/S; xq = Qo/X; and £ = E/S. The dimensionless potential strength Uq is 
related to the wall size iV = A w X/Iq by 

\n 2 g 2 fi 2 B J 

where g is the electron g-factor; and \Lb is the Bohr magneton. 



We utilize a shooting algorithm [113] to solve the eigenvalue problem numerically. Details 



of this calculation will be given in [[UJ. Our results for the tunneling gap A appear in 
Tables 1 and 2 and correspond to macroscopic and mesoscopic walls of thickness A = 1000A 
and 50A respectively. So far we have assumed our system of wall and voids to be completely 
isolated. If a weak stray magnetic field tt e xt were present, it would produce a bias e in the 



gap A = J Aq + e 2 . This bias is a consequence of the Zeeman energy density —~M-H ext . For 
an actual stray field, the direction of Ti ext is unknown and the experimentally relevant bias 
e is obtained by averaging over this direction. It is easily shown that e ~ MA w LH ext . We 
take H ext ~ 10 _6 G as indicative of the magnitude of a stray magnetic field (AC-magnetic 
fields of 10~ 5 G have been used in measurements of the frequency dependent magnetic 
susceptibility 0). Clearly, a necessary condition for observable domain wall QC is A > e. 
For N > 10 4 , we find that A w = (1.1 x l(T 17 cm 2 )iV (and M D = (5.6 x l(T 28 gm)iV), so 
that for LaGaYIG for which M ~ 10G, e = (1.04 x 10 _11 if)JV. Thus the bias grows 
with wall size N, making observation of QC more difficult for the larger walls. To proceed 
further, note that there exists a limited range of void separations for which Ao corresponds 
to QC and still satisfies A > e. For L < L min , the ground state energy is above the 
central barrier; while if L > L max , A < e. Let 1Z = L max — L min be the size of the 
allowed range of void separations. Its value is given in Table 1 and is obtained (for a 
given N) by comparing A and L with e. Quantum coherence will be observable only 
when: (i) the uncertainty in L satisfies AL <^ TZ; and (ii) 1Z ^> C, where C ~ (2 — 3)/ 
is the coarse graining length scale. One expects that AL ~ Iq and in our calculation 
Iq = 5A. If either (or both) of these conditions is (are) not satisfied one would expect that 
the experimentally relevant gap would be an average of Ao over the appropriate length 
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scale. From Table 1 we see that, for N > 10 4 , such an average is necessary and that any 
reasonable procedure gives A < e. Thus macroscopic QC (JV > 10 4 ) is not expected to 
be observable due to the rapid variation of the tunneling gap to small changes in L and 
the large bias e introduced by a stray magnetic field. We also see that in the case of the 
SDWP (in the absence of dissipation), the conditions for observable QC do not rule out 
walls with N ~ 10 2 — 10 3 (see Table 2) which would correspond to mesoscopic quantum 
coherence. (Here A w = (2.67 x l(T 16 cm 2 )iV; M D = (2.73 x l(r 25 gm)iV; and e is given m 
Table 2.) We now go on to examine the effects of asymmetry on the case of mesoscopic 
QC. 

Since U = 2KVd, the larger the void, the more strongly it attracts the wall. Thus, if 
asymmetry is sufficiently pronounced, QC is lost because the larger void pins the wall. This 
effect can be seen by following the ground state energy as we increase a from 1 (see Table 
3). Imagine a = 1 (corresponding to a SDWP) and that {N,U ,L} are such that, in the 
absence of dissipation, we have QC. Imagine further that we increase a so that the void at 
q = Qo attracts the wall more strongly than the other void. This stronger attraction causes 
E to decrease (i. e. become more negative) as the probability distribution in the ground 
state begins to shift towards q = Qo. As we continue to increase a, we reach a critical value 
a* at which E is equal to the value of the AsDWP at the metastable minimum of the left 
well U meta . For a > a*, E drops below the metastable minimum which corresponds to 
the pinning of the wall at q = Q and the destruction of QC. Intuitively, we expect that 
when AUo = Uq — U$ = (a — l)U ~ is approximately equal to the barrier height Ubh of the 
SDWP, the larger defect will pin the wall. For mesoscopic walls, U^h is given in Table 2. 
For N = 300, L = 75A; U bh = 2.8 x 10~ 7 eV (note the difference in units relative to Table 
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2). Thus AU = Ubh corresponds to aZ 1.008. A numerical calculation of E for this 
case gives a* = 1.038 > (see Table 3). For the spherical voids we have been considering, 
if R_ = 10 A, then R + = a* 1//3 i?_ = 10.1 A. Thus if asymmetry is not to destroy QC, 
the radii of the two voids must satisfy AR = R + — R_ < 0.1 A. Such a tolerance is 
clearly unattainable and since AR <Cwe must average the effects of asymmetry over the 
coarse-graining length scale C ~ 10 — 15A As the majority of AR values entering into the 
average correspond to pinning of the wall, we conclude that asymmetry in the voids acts 
to destroy QC in this case. We might hope to overcome this difficulty by increasing L and 
so increasing U bh . For L = 103A, N = 300, the SDWP tunneling gap is A = 7.7 x 10~ 9 K 
(see Table 2). At this separation, QC is marginally observable in the absence of dissipation. 
In this case U bh = 2.8 x 10~ 6 eV (see Table 2). Then AU = U bh gives a; = 1.085. We 
did not determine a* numerically for this case. In the previous example we saw that 
(a* — 1) ~ 5(a7— 1) so that we will estimate (a* — 1) ~ 10(a7— 1) for this case. This gives 
a* ~ 1.85. For R_ = 10A, AR ~ 2.3A. Again AR < C so that an average of the effects of 
asymmetry over C ~ 10 — 15A is necessary. As in the previous case, the majority of the 
AR values correspond to pinning of the wall so that we again conclude that asymmetry 
in the voids will act to destroy QC in this already marginal case. Larger values of L lead 
to A < e. We see that asymmetry of the voids will be sufficient to destroy any remaining 
vestige of domain wall QC — even in the absence of dissipation. The basic difficulty is that 
maximizing the tunneling gap requires a very small central barrier; so small in fact, that 
the most minute asymmetry in the two voids produces a bias in the double well which is of 
order of the height of the central barrier and so capable of pinning the wall (viz. destroying 
QC). We suspect that this is generally true of macroscopic QC: large objects require small 
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barriers which are easily removed by small imperfections in the experimental set-up. 

In this paper we have carried out a numerical analysis of domain wall quantum coher- 
ence in a uniaxial magnetic insulator with quality factor Q ^> 1 at T = 0. We find that 
QC on any scale larger than microscopic appears unlikely due to the combined effects of 
stray magnetic fields and asymmetry in the voids which are responsible for producing the 
double well potential seen by the domain wall. Our calculation assumed a flat wall al- 
though curved walls are expected in the thin film scenario. For this scenario, and for voids 
of given size, tunneling will only occur if L < L crit when curvature effects are included 
||. This is because the curvature energy acts to raise the minima of the DWP relative to 
the top of the central barrier, thus reducing U^. When L = L crit , the central barrier has 
disappeared and we are no longer in the QC regime. If L crit > L max , the tunneling gap 
becomes unobservable before curvature effects become significant. Otherwise, L cr a < L max 
and curvature effects act to reduce the range of void separation 1Z — > lZ new = L crit — L min 
corresponding to QC. For macro-walls, TZ was already small enough to rule out macro- 
QC so that decreasing 1Z — > 7Z n ew acts to strengthen this conclusion. For meso-walls, if 
7Z new < C, then curvature effects have reduced the range of QC sufficiently that stray mag- 
netic fields are expected to make meso-QC unobservable. Finally, if 7Z new > C, asymmetry 
in the void sizes can more easily destroy meso-QC since curvature effects act to reduce 
Ubh- Thus curvature effects are not expected to modify our conclusion that observation of 
meso- or macro-QC of domain walls appears unlikely. Our calculation ignores dissipation; 
a proper inclusion of its effects will also act to strengthen our conclusion. Although we 
have considered a particular type of defect, we expect our conclusion to also apply when: 
(1) the defects are such that the pinning potential is primarily due to a reduction in the 

9 



exchange and anisotropy energies of the wall (when pinned to the defect); and (2) the wall 
configuration M(y) depends only on the wall normal coordinate y. Under these conditions 
we continue to expect U(q) ~ (<//A 2 + K)V d sech 2 (q/ \) and our results to be indicative of 
the scale of QC to be expected. It should be noted that the phase coherences necessary 
for establishing quantum coherence are much more delicate than those necessary for es- 
tablishing quantum tunneling so that our work does not preclude a priori the possibility of 
observing macroscopic quantum tunneling of domain walls in magnetic insulators. I would 
like to thank Philip Stamp for useful discussions, T. Howell III for support, and NSERC 
of Canada for financial support. 
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List of Table Captions 



1. Tunneling gap A ; width 1Z of range of observable void separations; central barrier 
height Ubh] and bias e for a macroscopic Bloch wall moving in a SDWP. 

2. Tunneling gap A ; width 1Z of range of observable void separations; central barrier 
height Ubh] and bias e for a mesoscopic Bloch wall moving in a SDWP. 

3. Asymmetry parameter a; ground state energy E ; metastable minimum U meta of 
AsDWP; and barrier height Ubh for a mesoscopic Bloch wall with iV = 300 and void 
separation L = 75A. 
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N(spins) L(A) A Q (K) 



K(A) U bh (K) e(K) 



10 4 1320 AB ~ 6 1.04 x 10" 7 





1322 


4.0 x 10" 4 


1.8 x 10~ 2 






1324 


2.0 x 10~ 5 


3.6 x lO- 2 






1326 


5.0 x 10~ 7 


5.9 x 10~ 2 




10 5 


1318 


AB 


~ 2 


1.04 x 10~ 6 




1319 


1.8 x 10~ 4 


3.0 x 10" 3 






1320 


9.6 x 10~ 6 


6.7 x 10" 3 






1321 


2.0 x 10~ 7 


1.2 x 10~ 2 




10 6 


1317 


AB 


~ 1 


1.04 x 10" 5 




1318 


2.4 x 10~ 5 


7.9 x 10~ 4 






1319 


9.1 x 10~ 9 


3.0 x 10~ 3 






1320 


< 7.2 x 10- 11 


6.7 x 10- 3 




10 7 


1317.2 


AB 


< 0.1 


1.04 x 10~ 4 




1317.3 


2.0 x 10~ 5 


8.5 x 10~ 5 






1317.4 


8.2 x 10~ 6 


1.4 x lO" 4 






1317.6 


4.0 x 10~ 7 


3.0 x 10" 4 






1318.0 


2.2 x 10- 10 


7.9 x 10~ 4 





*AB = ground state Above Barrier 



Table 1: 
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N(spins) L(A) A (K) K(A) U bh (K) e(K) 



300 74 AB* ~ 30 5.2 x 10~ 9 

75 1.1 x 10~ 3 3.3 x 10~ 3 

80 2.5 x 10" 4 7.2 x 10~ 3 

103 7.7 x 10~ 9 3.2 x 10" 2 

110 2.4 x 10~ 10 4.0 x 10~ 2 

120 1.8 x 10~ 12 5.2 x 10~ 2 



3000 69 AB ~ 10 4.3 x 10- 

75 3.1 x 10~ 6 3.3 x 10~ 3 

80 6.1 x 10~ 9 7.2 x 10~ 3 



*AB = ground state Above Barrier 



Table 2: 
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a E (K) 



U m eta {K) 



U bh (K) 



1.000 -1.4638 x 10- 1 -1.4859 x 10" 1 3.33 x 10" 3 



1.010 -1.4719 x 10" 1 -1.4894 x 10" 1 2.96 x 10" 3 



1.020 -1.4814 x 10- 1 -1.4930 x 10" 1 2.59 x 10" 3 



1.030 -1.4914 x 10- 1 -1.4966 x 10" 1 2.23 x 10" 3 



1.038 -1.4997 x 10- 1 -1.4995 x 10" 1 1.94 x 10" 3 



Table 3: 
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